Investigation of Some Physical Properties of Accretion Induced 
Collapse in Producing Millisecond Pulsars 



Of 



Ali Taani 1 and ChengMin Zhang 1 

National Astronomical Observatories, Chinese Academy of Sciences, Beijing 100012, China 

alitaaniObao .ac.cn 

Mashhoor Al-Wardat 2 
Department of Physics, Al-Hussein Bin Talal University, P.O.Box 20, 71111, Ma'an, Jordan 

and 

YongHeng Zhao 1 

National Astronomical Observatories, Chinese Academy of Sciences, Beijing 100012, China 



ABSTRACT 

We investigate some physical characteristics of Millisecond Pulsar (MSP) such as magnetic 
fields, spin periods and masses, that are produced by Accretion Induced Collapse (AIC) of an 
accreting white dwarf (WD) in stellar binary systems. We also investigate the changes of these 
characteristics during the mass-transfer phase of the system in its way to become a MSP. Our 
approach allows us to follow the changes in magnetic fields and spin periods during the conver- 
sion of WDs to MSPs via AIC process. We focus our attention mainly on the massive binary 
WDs (M > I.OMq) forming cataclysmic variables, that could potentially evolve to reach Chan- 
drasekhar limit, thereafter they collapse and become MSPs. Knowledge about these parameters 
might be useful for further modeling of the observed features of AIC. 

Subject headings: Stars: neutron stars, white dwarfs, cataclysmic variables, x-ray binaries, fundamental 
parameters 



1. Introduction 

Cataclysmic variables (CVs) are short-period 
binary systems consisting of massive white dwarfs 
(WDs) primary that is accreting material via 
Roche-lobe overflow from l ow-mass secondary 



stars ( s ee i. e. Warner J 19951) 



(2002J); iMiddleton et all (|201lh ). so Accretion In 



Warner and Woudt 



duced Collapse (AIC) i s favo rable in these sys- 
tems ( van den Heuvel ( 2004) ). However, the 
AIC scenario has been proposed as an alterna- 
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tive source of recycled pulsars sufficient to ob- 

viate the difficulties with the standard model 

Bhattacharva and van den Heuvell ( 1991 ): Hurley et al 



2010h : iTaani et all (^Olj vT While other au- 
thors questioned the viability of an AIC origin 
for the recycled pulsars on the o retical and sta- 
tistica l g rounds (IRomanil (1 1990); van den Heuvel 
(|l995t )). IFerrario and Wickramasinghel (|2008a Bbh 
argued that the AIC channel can form binary 
Millisecond Pulsars (MSPs) of all observed types 
with B ~ 10 8 ~ 9 G, P ~ 20 ms and e - 0.44. 



Bhattacharva and van den Heuvell (|l99lh 



Lorimerl (|2008h : IWang et all (|201lh ) The 



re- 



quired conditions for the formation of an AIC in 
the case of a steady accretion were summarized in 
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Nomotol (|1984 Il987t ). Ilsern and Hernanj (Il994l) 
and in two recent works IChen et al.l (|201ll ) and 
Taani et al 1 (l2012h . However, unless the accretion 



rates are very high (> lO _7 M yr _1 ), the accreted 
(Hydrogen) matter will not go into steady nuclear 
burning on the surface of the WD, and hence its 
mass will never grow. In the situation of lower 
accretion rates, the accreted matter will burn in 
Super Novae (SN) type la, such that the WD will 
never reach Chandrasekhar limit (1.44 M ). 

Nowadays, there is an increasing amount of 
evidence which points to the fact that at least 
a fraction of low-mass X-ray binaries (LMXBs) 
should be the r esult of AIC of a WD in a low-mas s 
bin ary sy stem (jTaam and van den Heuvell (Il986lh 
Li (2010)). However, this straightforward hypoth- 
esis had to face several problems, some of which 
are not yet completely solved. The behavior of the 
external layers of an accreting WD depends basi- 
cally on the rate of mass transfer and on the chem- 
ical composition of the accreted material. Even 
though a lot of work has been done on this issue, 
much of it is based on very restrictive assumptions 
(steady spherical accretion, mass for instance). 

In general, accretion of material lighter than 
carbon will induce its burning. Only steady burn- 
ing or non-disruptive flashing r egimes can be al- 
lowed in NS formation scenarios (jWickramasinghe and 



mechanism is similar to LMXB. While long or- 
bital periods could potentially occur via explosion 
power from supernova type la, which, circum- 
stances has no observational counterparts, con- 
sequently, create single MSP. It is currently un- 
known whether such AIC process appears and how 
much it contributes to form MSPs population, this 
is a focus of much present-day research. 

The purpose of this paper is to demonstrate 
how to infer some of the observable quantities and 
relate them to the dipolar magnetic magnetic field 
at the surface of the NS, bottom field strength 
(when the magnetospheric radius Rm = R) and 
spin period, during the AIC of a WD in its way to 
become a member of MSP family. 

To accomplish this goal, the model developed 
bv lZhang et ail (j2009h is extended to include the 
magnetic fields and radii for cataclysmic variables 
(CVs) observed sample. Thus a significant reduc- 
tion of the polar field strength could occur if such 
a system accumulates the required critical mass 
M cr u. Thus we demonstrate the observed prop- 
erties of the magnetic CVs, that can generally be 
explained by the model where the field is at most 
only partially restructured due to accretion. 



Ferrariol 



2005) . For which conditions this actually happens 



is still largely undetermined, but it seems that 
there are several combinations of the parameters 
mentioned above which allow the star to retain 
the accreted material and to grow in mass. Even 
though the WD manages to accrete mass steadily, 
the subsequent collapse is not yet guaranteed 
because explosive thermonuclear burning always 
happens in the contracting core before reaching 
Chandrasekhar limit. The burning might prop- 
agate through the entire star, releasing enough 



energ y to blow it apart (jWickramasinghe et al 
2009h . 




Observable parameters of binary MSPs, e.g. 
mass of the pulsar, mass of the companion, spin 
period of the pulsar, orbital period, orbital eccen- 
tricity, etc., are used to probe the past accretion 
history of the MSPs. 

AIC produces WDs having mass less than 
Chandrasekhar limit, because at least ~ 10% of 
the mass accreted by the W D goes to the b inding 
energy of the nascent NS ( Bagchil ( 2011 )); this 
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Fig. 1. — Mass distribution of all CVs observed 
sample. The solid line is the curve fitted using a 
Gaussian function. The data are taken from Ritter 
& Kolb (2011) to construct this distribution. 
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2. Methodology 
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Fig. 2. — Histograms of the bottom magnetic field 
in CVs. The distribution is relatively Gaussian. 
Here the field strength of most CVs has been in- 
ferred on assumption that CVs of non-spherical 
accretion on to magnetized NS wit h a centered 
dipolar field (see IZhang et al.l <|2009h ). The solid 
line is the curve fitted using a Gaussian function. 



2.1. The Observed Population of CVs 

Most of the behavior of CVs can be understood 
in the context of a close binary with mass trans- 
fer from a low-mass, near-main-sequence compan- 
ion star. The mass donor is known as the sec- 
ondary, and the star which accretes mass is known 
as the primary. In CVs the primary is a WD, 
and the secondary is usually a late-type roughly 
main-sequence star. Interacting binaries are inter- 
esting and important to observe, as they allow the 
study of accretion on observable timescales. More- 
over, much can be learnt about binary evolution 
from them, in particular the effects of accretion 
on their ultimate evolution. It is predicted that 
above the 2 — 3 hr orbital period gap, angular 
mo mentum losses are driven by magnetic brak- 
ing ( Ivanova and Taam (2003))'. while below the 
gap the systems evolve via the emission of gravi- 
tation a l wave radiation (see i.e. ISpruit and Ritter 
(|l983l ): iKnigge et ail (|201lh l. However, the theo- 
retical relationship between the mass and radius 
of a WD is essential in interpreting some observa- 
tional results. M-R relation was first defined by 
Chandrasekharl B as: 



Derived Observed 




Log B (G) 

Fig. 3. — Observed magnetic field distributions 
of MSPs (shaded histogram) and derived ones of 
the CVs (solid histogram), which were calculated 
by assuming the conservation of the magnetic flux 
(see the text). The da ta for MSPs were taken from 
the ATNF catalogue (jManchester et al.ll2005h . 
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f 



MVS 



(1) 



for Sirius B ( M = 0.986M Q , R = O.OO846R ), 
where mass and radius are measured in standard 



solar radii and masses ([Woodl (|1995[ )). We use this 
M-R relation to derive the radii of CVs which are 
not yet observed. The results of the calculated 
M-R values for massive CVs are listed in Table Q] 
The sample of mass of CVs we have consid- 
ered is the set of bin ary systems collected by 
Ritter and Kolbl ( 201ll ) from several high-quality 
surveys. The sample consists then of 133 systems. 
Fig. 1 shows the Gaussian distribution of mass of 
CVs, with mean at M C v ~ 0.822M Q , where the 
median mass of these sources is Mcv = 0.81 ± 0.16Mq. 
Since the Ritter and Kolb (2011) catalog gives 
a mass of CH UMa M cv ~ 1-9 ± 0.30M© (well 
above the Chandraskhar mass limit for a WD), we 
decided to exclude this source from our study. We 
fit the Gauss function to the mass distributions. 
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The Gauss function we choose reads, 

y = y 0+ *ew(-2(^) 2 ) (2) 



i/'o 



where yo, x , wq and A are offset of y-axis, center 
of x-axis, width and area represented by the curve 
respectively. The fitting results are listed in Table 
2. 

2.2. Magnetic field 

The strength of the magnetic field influences 
the physical conditions in the accretion column, 
and as a consequence the dominate radiation 
mechanisms and the e volutionary phase of CV s 
population (see i.e. 



Spruit and Rittc 



te3 
)83): 



1983); 



Qian et al 



King et all (|1994l ): lRappaport et al.l(| 198 
(|201lh ). To derive the bottom magnetic field 
in CVs (where the field is at most only par- 
tially restru c tured due to accretion), we follow 
Zhang et all (|2009h model, which suggest that the 
bottom filed Bf for WDs is reached when the 
magnetospheric radius Rm = R, giving 

B f WD = 2.8 x 10 3 M 16 1/2 m 1 / 4 i?~ 5/4 <T 7/4 (G) 

(3) 

where, m — M/Mq is the WD's mass in CVs 
M, where the mean mass of these sources is 
Mcv = O.822M0, R 9 is the WD radius in units 
of 10 9 cm, <f> is a parameter which is estimated 
to be ~ 0.5 for disc accretion, and M±e is the ac- 
cretion rate in units of 10 16 gs _1 . The concept is 
displayed pictorially in Fig. [2j which shows clearly 
two peaks (bi-modal) Gaussian distributions for 
most of the CVs with magnetic fields in the range 
(0.60 - 11.20) xlO 3 G. We fit the Gauss function 
(same as Eq. (2)) to this distribution. A full sum- 
mary of the main fitting results is given in Table 
3. 

Another interesting point is to investigate the 
correlation between the magnetic fields of MSPs 
with bottom fields of CVs, which can be derived 
according to the conservation of the magnetic flux, 
in other words, study the magnetic fields that pro- 
duced by AIC. Let's start, then, with the conser- 
vation of the magnetic flux 



B 



NS 



Bf.WD x 



/ Rwd \ ' 
V Rns J 



(4) 



As an example of a typical NS we suggest to 
adopt a radius (Rns) of about 15 x 10 5 cm, where 
a can be expresses as (Rwd /Rns), subsequently 
a 2 ~ 10 5 ' 5 ~ 6 , and B f ~ 10 3 G in CVs. As a result, 
one arrives at an expression for the minimum value 
of magnetic field in NSs equals to: 



B 



NS 



10 8.5-9 G 



(5) 



which is somehow similar to the observational val- 
ues of MSPs. 

To study the contribution of magnetic field of 
CVs on the MSPs population, we produced a his- 
togram of the total sample of MSPs (Fig. 3), 
which shows a clear Gaussian distribution for most 
of the objects, while the observed sample shows 
two peaks and a small fraction of these systems 
have intersection. We note that the pinnacle of de- 
rived histogram is much higher than in observed, 
this is due to the sensitivity of the radio obser- 
vation process for MSPs. Since all pulsar surve ys 
have some limiting flux density ( Lorimer ( 20081 )). 
However, new data acquisition systems provide 
significantly improved sensitivity to MSPs and are 
being used to resurvey the sky. Therefore the ob- 
served numbers of MSPs of magnetic field are rel- 
atively high compare with CVs population. It is 
noteworthy to mention here that we also found 
that ~ 15% of MSPs could potentially be formed 
from CVs via AIC process. The peak of the dis- 
tribution for derived sample occurs at ~ x 10 9 4 G 
comparing with ~ x 10 8 5 G for the observed one. 
Furthermore, very high magnetic field Polar sys- 
tems (like AM Her) are not represented in this 
figure. 

2.3. Spin period 

It is well established that the NS in binary 
MSP systems forms first, descending from the ini- 
tially more massive of the two binary stellar com- 
ponents. The NS is subsequently spun- up to a 
high spin frequency via accretion of mass and an- 
gular momentum once the secondary star evolves 
(For a more detailed review, see i.e. Alpar et al.l 
( 1982);|Bhattacharva and van den Heuve ll (Il99ll): 
Tauris et al.l (|2000l ): iLiu and Chenl (|201ll ): iTaurisl 
(|201ll )) In this section we will determine the spin 
period of MSPs originated from WDs, the process 
begins with simple Keplerian frequency, the an- 
gular velocity of the NS is equal to the Keplerian 
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angular velocity of the magnetosphere, at roughly 
the Alfven surface, 



lj-3/2 



P. 



MSP 



R 



3/2 
NS 



(6) 



from which we obtain Pmsp a s a function of Pwd 



MSP 



W D ,min 



R 



NS 



R 



WD 



3/2 



(7) 



where t he minimum spin pe r iod of the stan 



dard W D (|de Jager et all (|1994[ ): lMiddleton et al 



(2011)) taking value ~ 30 s, Rns — 10 km and 
Rwd — 1000 km (see Eqn 1) , giving 



P 



MSP 



lms 



Note that the estimated spin periods are not 
far from the true spin periods of the MSPs. This 
clearly supports the view that the AIC effect play 
an important role in deriving spin periods for 
MSPs. 

Fig. 4 shows the observed and derived spin pe- 
riod distributions of MSPs. It can be found that 
most of derived values centered at 0.011 s and 
ranging from 0.011 to 52.8 s. As for the observed 
MSPs, shows relatively Gaussian and regular dis- 
tributions for most objects, centered at 4.4 ms and 
ranging from 1.4 ms to 19.5 s. According to these 
distributions, the ratio of MSPs originated from 
CVs is about ~ 10%. This result agrees wi th some 
theore tical predictions such as those by IWarner 
(|l995h and lWarner and Woudtl (|2002h . 



2.4. Mass 

Most of the stellar binary systems consist of 
two well-separated components, with a little or no 
interaction beyond their mutual gravitational at- 
traction. However, in close binary systems tidal 
effects become important, leading to a circulariza- 
tion of the orbit, and spinning up the star's ro- 
tational speeds so that they become synchronous 
with the orbital motion. 

The current edition ( 7.16)Q] of the Ca t alogue 
of Cataclysmic Binaries iRitter and Kolbl (|201ll ) 
contains a wide variety and natures of sources 
such as 880 CVs, 98 LMXBs, and 319 related 
objects. Among them we have just 18 mas- 
sive CVs in the range 1.1 — 1.3M Q , which is 



1 http: / /www. mpa-garching.mpg.de/RKcat/ 



n 



R 



] Derived 
] Observed 



R 



1 1 



V 



v 



-2.8 -2.6 



L °9 P Spln [S] 



Fig. 4— The distribution of MSPs on the 
base of spin-period; Observed (shaded histogram) 
and derived (solid histogram). The data for 
MSPs were taken fro m the ATNF catalogue 
([Manchester et al.ll2005h 



considered a relatively small number. A sum- 
mary of the known properties of these sys- 
tems is given in table [TJ. If matter is ac- 
creted at a rate of M ~ 10~ 9 M W /yr to 10" 8 M Q /yr 
( Zhang and Koiima (2006)), and the total mass 
accreted exceeds a critical value AM crit ~ 0.1 — 
0.2M Q with 10 9 yrs, then the massive CV will 
be recycled to become a MSP after the mass 
reaches the Chandrasekhar limit. It is worth 
noting here that the mass of the newly formed 
object will be be l ow th e maximum mass for MSP 
van den Heuvell (|2004h ). This provides evidence 



for the AIC scenario in massive CVs and evolu- 
tionary hypotheses of MSP birthrate. In sum up, 
whether a CV might possibly do this will depends 
on: (i) the mass of the WD in the system, and 
(ii) the mass transfer rate that the donor star may 
supply, because we know that at low mass transfer 
rates one only gets SN explosions Type la, and the 
WD will never grow to become a NS. 

It should be noticed that the AIC process leads 
to MSP with mass less than Chandrasekhar limit. 
When a WD gains mass by accreting matter from 
its companion, its binding energy also increases, 
this effect is significant in the estimation of the 
amount of mass accretion. However, unfortu- 
nately most of the studies argued that the bind- 
ing energy of the NS is not commonly considered 
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Table 1: List of physical properties for the all massive CVs (M > 1.0Af Q ) * 



Name 


M/Mq 


Porb (D) 


RDerivedi^Cm) 


BDerievd (lO^G) 


CAL 83 


1.3 


1.04 


0.76 


13.95 


CI Aql 


1.2 


0.618 


0.83 


12.38 


BV Cen 


1.24 


0.61 


0.8 


13 


EY Cyg 


1.1 


0.49 


0.9 


10.86 


RU Peg 


1.21 


0.37 


0.82 


12.53 


QZ Aur 


1.05 


0.357 


0.95 


10.13 


MU Cen 


1.2 


0.342 


0.83 


12.38 


BF Eri 


1.28 


0.27 


0.78 


13.63 


VY Scl 


1.22 


0.232 


0.81 


12.69 


RX And 


1.14 


0.209 


0.87 


11.46 


SS Aur 


1.08 


0.182 


0.92 


10.57 


BD Pav 


1.15 


0.179 


0.86 


11.61 


U Gem 


1.2 


0.176 


0.83 


12.38 


V603 Aql 


1.2 


0.138 


0.83 


12.38 


WW Hor 


1.1 


0.08 


0.9 


10.86 


CU Vel 


1.23 


0.078 


0.81 


12.84 


DP Leo 


1.2 


0.062 


0.83 


12.38 



*The masses and orbital periods were taken from lRitter and KolrJ d201lT) . while B and R are determined by the condition that the 
magnetospherc equal to the WD radius (Zhang et al. 2009). The model-dependent parameters are <p = 0.5 and M = 10 16 gs — 1 , 
(see the text). 



Table 2: The fitting results of distributions for mass of Cvs. 

quantity tjq ~x~ c w A ~o x^/Dof R 2 



M CVs 0.53 ±0.64 0.79 ±0.012 0.40 ± 0.03 12.30 ± 0.92 0.20 3.311 0.955 



Table 3: The fitting results of distributions for derived magnetic field strength in CVs . 



quantity 


yo 


X Xi Wo w 


A 


B R l 


^Derived 


-0.05 ±0.004 


3.93 ±0.04 3.41 ±0.08 0.186 ±0.008 0.37 ±0.003 


14.42 ±4.57 


5.45 ± 1.35 0.942 
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(|Bagchil (120111) 1 



2.5. Orbital periods and eccentricity 

We now believe that in the case of spherical 
SN explosions in which C-0 WD (at low mass 
transfer rates) produced impart significant kicks to 
the companion that make their survival prospects 
within binary systems in a long orbital period 
or bleak. The asymmetric mass loss during the 
AIC process is expected to provide kick velocities 
of abo ut 50 kms -1 or less (jGothoskar and Gupta 
(2000])). Consequently, the AIC produces very low 
eccentricities (e < 0.1). When the natal MSPs re- 
ceive a relatively strong kick (> lOOkms -1 ), the 
AIC channel produces eccentric binary MSPs in 
the Galaxy. Such a kick seems to be highly un- 
likely in the AIC process, hence the probability of 
forming eccentric binary MSPs via the AIC chan- 
nel can be ruled out. Even if a high kick is allowed, 
the AIC channel cannot produce ecc e ntric with an 
orbital period of > 70 d lchen et all (|201ll ). 

In sum up, the high-mass binary systems tend 
to have highly eccentric orbits, while the low-mass 
binaries have smaller eccentricities. The origin of 
MSPs in low-mass binaries could be due to either 
AIC or supernova explosion. If MSPs are formed 
by the AIC process, their average velocity would 
be much lower than that of normal pulsars. Fur- 
ther work will certainly clarify this issue as known 
sources are better characterized. 

3. Summary and Conclusions 

We have studied the possible AIC scenario in 
massive CVs, through simple calculations of some 
physical characteristics such as magnetic field, 
spin period and mass. The basic information that 
goes into the calculation can be summarized as 
follows 

1. CVs would be invoked via the capability 
of producing a significant portion of the 
MSPs population through the AIC process, 
a regime which otherwise may be unattain- 
able by the normal channels that produce 
NSs (i. e. long orbital periods for half of the 
MSPs and the formation of single MSPs) in 
order to reduce the gap between standard 
LMXB evolution and the observed popula- 
tion of MSPs. 



2. We find that the quantitative implications 
of our calculations are that, we can estimate 
the expected magnetic field strength in the 
observed population of MSPs which could 
be contributed from CVs is about ~ 15%. 
While it turns out that the contribution 
in spin period is about ~ 10%. It would 
be very interesting though, if observations 
would yield either a very slowly spinning 
(few hundred ms) pulsar associated with a 
very low magnetic field, or a high magnetic 
field MSP with a very rapid spin. (However, 
the latter kind of pulsars might be unlikely 
to be detected given that its very strong 
magnetic dipole radiation would slow down 
its spin rate within a few Myrs). 

3. We compute the bottom magnetic field in 
CV samples, with a dis cussion on th e impl i- 



(|2009f ). 



cations of the model in IZhang et al. 
We extend this model to conditions appro- 
priate to accreting WDs. Thus the model 
predicts that the magnetic polar cap widens 
as material is accreted with the field being 
advected towards the equatorial regions by 
the ensuing hydrodynamic flow. 

4. On the basis of our results we conclude that, 
it is possible to make comparisons about the 
relative importance of spin periods, mag- 
netic field and some other parameters in CV 
population, that potentially produced via 
AIC process about their aggregate impact 
in the observed population of MSPs. 

5. The mass of MSPs which produced via AIC 
process is less than Chandrasekhar limit, be- 
cause at least ~ 10% of the original mass 
goes to the binding energy of the NS. 

6. We further find that the predictions of some 
parameters after AIC process for the aver- 
age levels are consistent with the observed 
population (such as spin period, magnetic 
field and mass). Future work will go steps 
further, and consider other observable quan- 
tities such as orbital period, eccentricity and 
mass ratio (q) for MSP populations using the 
more extensive data se t now available from 
Ritter and Kolbl (120111) and other sources. 



Ultimately, the critically unique analysis will 
be the cross correlation with the CVs, AIC pro- 
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cess and MSPs, which is well beyond the scope of 
this letter. Our main conclusion is that the simple 
calculations presented herein, based on the stan- 
dard relationships between mass, magnetic field, 
spin and orbital periods, required to support the 
more sophisticated treatments of the problem, and 
suggests that how the AIC contributes to the ob- 
served population of MSPs. 
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